Suppression of Complex Spiral— Wave Activity in an Ionic Model of Cardiac Tissue by 

Weak Local Stimulations 
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On the basis of a quite realistic ionic Fenton-Karma model of the cardiac tissue we consider the 
problem of defibrillation by a local weak forcing. In contrast to other systems, this model accu- 
rately reproduces the most essential mesoscopic properties of the cardiac activity. It is shown that 
suppression of spiral-wave turbulent dynamics in the heart tissue may be realized by a low-voltage 
local non-feedback electrical stimulation of monophasic and biphasic shapes. After stabilization the 
medium goes to a spatially homogeneous steady state. 

PACS numbers: 05.45.Gg, 47.54.+r, 82.40.Bj 



Cardiovascular diseases (CVDs) are responsible for 
more than 4 million deaths each year in Europe (over 
1.5 million deaths in the EU) and account for about 30% 
of life-span loss in Europe (over 30% in the EU) 0. A 
major group of CVDs involves disturbances of the nor- 
mal cardiac rhythm (cardiac arrhythmias). The extreme 
form of cardiac arrhythmias and the prevalent mode of 
the sudden death among patients with CVDs is ventric- 
ular fibrillation (VF), which is a fast developing distur- 
bance of spatially organized contraction of ventricles that 
is a consequence of abnormalities of electrical conduction 
in the heart muscle. Following to the contemporary con- 
jecture, VF is produced by a multiple wavelet re-entry, 
which is spiral waves in 2D and scroll waves in 3D Q (i.e. 
by spatio-temporal chaos or spiral-wave turbulence). 

Since VF, sustained for only a few minutes, leads to 
death, an immediate intervention is required. In emer- 
gency care medicine the application of high-energy elec- 
trical stimulation through the patient's chest is com- 
monly used to suppress the fibrillation and restore the 
normal rhythmicity of the heartbeat. However, high- 
energy shock can cause the necrosis of myocardium or 
give rise to functional damage manifested as disturbances 
in atrioventricular conduction. 

The application of electrical pulses for the termination 
of fibrillations is also used in implantable cardioverter 
defibrillators (ICDs). These devices surgically implanted 
into the bodies of high-risk cardiac patients and initi- 
ating low-power electrical pacing pulses automatically 
when they detect a dangerous activity. However, in the 
case of complex arrhythmias (but not fibrillation) , some 
patients may also undergo the ICD action. As a con- 
sequence, such patients additionally have a severe pain. 
Therefore, a very important factor in the design of mod- 
ern ICDs is decreasing the stimulation amplitude in order 
to avoid a painful high energy shock and damage to the 
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heart itself and surrounding tissues. Thus, there is high 
demand in clinics on alternative methods of defibrillation 
which would work with lower voltages. 

The recent research Q may provide an alternative to 
the conventional ICD therapy by terminating re-entrant 
arrhythmias with the field strengths that are 5 to 10 times 
(or delivered energy 25-100 times) weaker than usual de- 
fibrillation shocks. However, the method used in |3j is 
valid only for the high-risk cardiac patients who had pre- 
vious myocardial infarctions ("heart attacks"). 

Theoretical studies suggest that low-energy defibrilla- 
tion protocols are also possible at exploiting the dynami- 
cal properties of re-entrant waves under electrical forcing, 
known as feedback-driven resonant drift Q. However, 
the major problem in practical use of the resonant drift 
is the change of the resonant frequency with the posi- 
tion of rotating wave, in particular, close to unexcitable 
boundaries. 

The qualitatively different approach of low-amplitude 
suppression the complex dynamics of nonlinear systems 
by application of perturbations without feedback was 
firstly proposed in |5j and mathematically substantiated 
in Recently this method was tested on a model of 
the Belousov-Zhabotinsky reaction and showed its valid- 
ity |3- The same parametric controlling but involving 
the feedback was used in Q . Although parametric (non- 
feedback) suppression or (feedback) controlling leads to 
the stabilization of complex dynamics, it is not realiz- 
able for the electrical defibrillation. The only application 
seems to be in formulating drug therapies which modify 
ionic currents in order, for example, to prevent alternans 
(beat-to-beat alternation in the action potential dura- 
tion) which is presumably one of the causes of breakup 
of a single rotat ing wave into multiple re-entrant waves 
(see, e.g. HEEf). 

The recent investigations of the active medium theory 
offer new opportunities for the electrical defibrillation: 
The amplitude of the external stimulation can be essen- 
tially decreased and the turbulent regime in excitable 
systems may be suppressed by a sufficiently weak pe- 
riodic external forcing applied globally [13 or locally 
[IH EH US [HI By this manner, it is possible not only to 
suppress spatio-temporal chaos and stabilize the media 
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dynamics, but also reestablish the initial cardiac rhythm, 
because after stabilization the medium goes to a spatially 
homogeneous steady state. 

In the present paper, on the basis of a realistic cardiac 
model we resolve the problem of suppressing the fibrilla- 
tive activity by a low-voltage local non-feedback electri- 
cal stimulation of monophasic and biphasic shapes. In 
contrast to FitzHugh-Nagumo type models, this model 
accurately reproduces such mesoscopic characteristics as 
action potential duration (APD), restitution and conduc- 
tion velocity (CV) found in cardiac tissue. It turns out 
that APD (CV) is a function of both the previous APDs 
(CVs) and the time between excitations, also known as 
the diastolic interval (DI) or recovery time 0] • 

In our investigations we used a three variable simplified 
ionic model (SIM) of the cardiac action potential, so- 
called Fenton-Karma equations 0, 0] : 

d t u = V(DVu) - {Jfi(u, v) + J so (u) + J sl (u, w)) , 
d t v = Q(u c —u)(l - v)/t~(u) - 9(w - u c )v/t+, 
d t w = 9{u c - u)(l - w)/t~ - Q(u - u c )w/t+, 

(1) 

where u is a dimensionless membrane potential; v,w are 
a fast and a slow ionic gates, respectively; D is a diffusion 
tensor which is a diagonal matrix in our case (isotropic 
simulations); J/j, J so , J S i are scaled ionic currents corre- 
sponding to the Na, K, Ca currents, respectively; Q(x) 
is a standard Heaviside step function (other functions, 
equations for the ionic currents and parameter values 
see in |3). The parameters correspond to the steep 
APD restitution (fitted to accurately represent the APD 
restitution in the full Bceler-Reuter model) with breakup 
close to the tip. 
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FIG. 1: The number of PSs as a function of time for the 
SIM with parameters corresponding to the Set3 in the original 
Fenton-Karma model 

The SIM supports many different mechanisms of spi- 
ral wave breakup into complex re-entrant activity. Our 
numerical simulations were performed in a 2D grid of 
500 x 500 elements corresponding to the tissue size of 
12.5 x 12.5 cm. We used periodic boundary conditions, 
which correspond to the torus topology. This geometry is 



more close to a real geometry of ventricles than a sheet 
of tissue (Neumann boundary conditions) and allowed 
us to exclude wave attenuation. As a measure of the 
suppression of turbulent effectiveness it is convenient to 
compute a number of phase singularities (PSs), i.e. tips 
of re-entrant waves, which can be detected by a num- 
ber of various techniques [T^. l20l | . We used the method 
described in |2*lj . 

During numerical analysis we have found that an ini- 
tial archimedean spiral wave breaks into complex tur- 
bulent pattern after approximately 500 ms (see Fig^). 
Chaotic system states of 500, 600, 700 ms were consid- 
ered as initial ones for all our suppression attempts. To 
hnd the suppression effect, we added the external almost 
point periodic forcing J ext (t) of the frequency and the 
amplitude A to the media. So, to stabilize the spatio- 
temporal chaotic dynamics we generated in the medium 
a single pacemaker (external electrode). 
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FIG. 2: Monophasic and biphasic impulses. 

The shape of external stimulation J ext (t) is one of the 
key factors strongly influencing on the suppression ef- 
fectiveness. The defibrillation shocks used in the clini- 
cal practise are of rectangular monophasic and biphasic 
shapes. In the contrast to defibrillation by pulses (a sin- 
gle shock or series of shocks), applied to the entire muscle 
or quite large part of it, we applied periodic stimulation 
of the same mono(bi)phasic waveforms (Fig|5J to a point 
of a medium (2x2 nodes). Optimal values of r (r < 1) 
vary as 0.05 -j- 0.15 for monophasic and 0.7 -r- 0.75 for 
biphasic stimuli. We took r = 0.1 for the monophasic 
stimulation and r = 0.7 for the biphasic waveform. 

The important problem is to select stimulation am- 
plitudes. To determine the highest possible values of 
A corresponding to experimentally observable shock- 
induced variations of membrane potential, we applied sin- 
gle stimuli of various amplitudes and durations. For the 
monophasic waveform duration of influence, Tt was 40 
ms, which corresponds to T — 400ms, r = 0.1 for the pe- 
riodic stimulation. For the biphasic stimulation duration 
of the positive influence, T(l— r) was 100 ms (T = 333ms, 
t = 0.7). FigEO shows membrane potential at a pace- 
maker site as a function of time for both stimulation 
shapes and three amplitudes. To ensure experimentally 
observed values of membrane potential for suppression we 
took A — 10^,A/cm 2 for both waveforms. Larger stim- 
ulation amplitudes produce the positive shock-induced 
variations of membrane potential larger than the maxi- 
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FIG. 3: Membrane potential at a pacemaker site vs time under excitation by single stimuli with amplitudes 10, 20, 30/iA/crn 2 
(bottom up). Left-hand side — monophasic stimulation, right-hand side — biphasic stimulation. 




FIG. 4: The frequency uiout of target waves as a function of the pacemaker frequency uim- Left-hand side — monophasic 
stimulation, right-hand side — biphasic stimulation. 
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FIG. 5: The number of PSs as a function of time for the 
SIM during monophasic stimulation of u)i n — 7 Hz and A — 
W/iA/crn 2 . Suppression onset are: t = 500 ms (solid lines), 
600 ms (dashes). 



mal observed value of 100 mV j^. Because the cardiac 
membrane potential varies as —90 -j- +30 mV, the maxi- 



mal positive value of membrane potential after applying 
shock is about 130 mV. Determination of the lowest am- 
plitude corresponding to the effective suppression is a 
task of the future research. 

We have measured the APD restitution curves ob- 
tained by two successive SI and S2 stimuli of various am- 
plitudes and durations (we use the 80% cutoff (APD S0 ) 
when calculating restitutions). It was found that the 
APD restitution curves are the same for all stimuli and 
exactly look like one in Fig. 4 in |l8|. However, the S1-S2 
interval of various stimuli is different to get the small- 
est DI of 43 ms. It is less when stimulation amplitudes 
and/or durations are greater. This means that the stim- 
ulation of larger amplitude and/or duration excites the 
heart tissue being in relative refractory state earlier. This 
inverse dependence of the refractory period on ampli- 
tude and duration of excitation phase does not appear in 
FitzHugh-Nagumo type models. 

Consider the problem of detecting the excitation fre- 
quencies, which provide the effective stabilization. Since 
a search of the suppression frequencies at random is in- 
effective, it is necessary to select the frequency of the 
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FIG. 6: Time evolution of the excitation pattern during biphasic stimulation started at 600 ms with close frequencies. (A) 
Effective suppression of complex activity by biphasic stimuli, u>in = 7.25 Hz. At the subsequent forcing all sites recovered 
without further activation. (B) Recovery of spiral-wave activity after its suppression by biphasic stimuli, u>in = 7.0 Hz. (C) 
Unsuccessful attempt to suppress turbulence by biphasic stimulation, Wj n = 6.75 Hz. 



stimulation close to the maximal possible frequency for 
a given medium. 

To find such frequencies, one can measure the period 
of target waves emitted by the created source (electrode) 
as a function of its own period and then choose values in 
the frequency intervals near the maxima of the frequency 
dependence. As is known, if there are several co-existing 
sources of periodic waves in an excitable medium, the 
interaction of waves leads to suppression of the sources 
with a longer period by a source with a shorter period . 
This was firstly formulated in |2^| and is caused by the 
destructive interaction of colliding waves in media, which 
mutually annihilate. If the leading source is an external 
electrode, it can suppress re-entry subject to the correct 
choice of its frequency and shape of stimulation. 

Thus, to find the suppression frequencies we generated 
pacemakers in quite small media volumes and determined 
the frequency uj out of the target waves as a function of the 
pacemaker frequency u>i n (FigQJ. One can expect that 
the spatio-temporal chaos can be suppressed by the point 
external perturbation in the frequency intervals corre- 
sponding to maxima of these dependencies. 

First, we tried to suppress complex activity by 
monophasic stimulation with u>i n = 3.13, 7, 10 Hz, cor- 
responding to the frequency maxima on the left-hand 
side of Fig^J It was found that although the suppres- 



sion phenomenon was observed for u>i n — 3.13 Hz and 
LUi n = 7 Hz, it strongly depends on the suppression on- 
set. For example, the monophasic forcing of u>i n = 7 Hz 
leads to the stabilization of chaotic dynamics if it starts 
at 500 ms, but suppression is unsuccessful if the suppres- 
sion onset is 600 ms (Fig|5} . Vice versa, stimulation with 
LUi n = 3.13 Hz started at 600 ms was successful. Suppres- 
sion of turbulence by stimuli with coi n — 10 Hz failed for 
all considered three initial conditions (500 ms, 600 ms, 
700 ms). This phenomenon is connected with the initial 
orientation of spiral tips with respect to the excitation 
source and remans to be explored. 

Second, we forced the system by the biphasic stimu- 
lation with uji n — 3.13, 7.25, 10 Hz corresponding to 
the frequency maxima on the right-hand side of FigQ] 
Again, the perturbation with iOi n — 10 Hz was unsuc- 
cessful. But, in contrast to the monophasic stimulation, 
biphasic forcing leads to the stabilization of complex dy- 
namics by stimuli with iOi n = 3.13 Hz started at 500 ms 
and by stimulation with u>i n — 7.25 Hz started at 600 ms. 
Because the frequency interval corresponding to the sec- 
ond maximum on the right-hand side of Fig0| is a quite 
wide, it is rather complicated to select the appropriate 
value of the stimulation frequency. Fig. [5] shows suscep- 
tibility to its choice. There is time evolution of the exci- 
tation pattern during biphasic stimulation started at 600 
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FIG. 7: The number of phase singularities as a function of 
time for the SIM during biphasic stimulation of LOi n = 7 Hz, 
A — 10/i A I err? (solid lines), 20fiA/cm 2 (dashes), 30/iA/cm 2 
(dots). The suppression onset is t = 600 ms. 

ms with close frequencies 7.25 Hz (panel A), 7 Hz (panel 
B) and 6.75 Hz (panel C). The panel A corresponds to 
the successful suppression (just an external pacemaker 
remained), B — to the recovered turbulence, C — to 
the unsuccessful suppression. As it was predicted, the 



patterns of the suppressed and recovered turbulence un- 
til suppressing spiral-wave activity were similar due to 
the small difference in the stimulus length, but creation 
of new spirals from an external pacemaker in the latter 
case (panel B) was unexpected. 

It should be noted that stimulation with u>i n = 7 Hz 
started at 600 ms (resulting in recovered turbulence when 
A = lQiiA/cm 2 ) leads to the effective suppression at dou- 
bling the stimulation amplitude. However, trebling the 
amplitude is not helpful (Fig Q). So, there is a nonlinear 
dependence on the stimulation amplitude. 

Thus, although the suppression effectiveness strongly 
depends on the stimulation frequency, amplitude and ini- 
tial conditions (suppression onset), it was found that for 
the correctly chosen these values the re-entrant waves can 
be easily eliminated. It should be noted that the ampli- 
tude of the excitation is by three orders of magnitude 
less than used in ICDs. Such low-voltage defibrillation 
has a great advantage because it does not require knowl- 
edge of the re-entry frequency. Moreover, in the case of 
VF, all the rotating waves are suppressed simultaneously 
and the initial cardiac rhythm can be reestablished. By 
these reaosns, probably the described new defibrillation 
strategy may be realizable in practice. 
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